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Abstract 

Using the entire OPAL LEP1 on-peak Z hadronic decay sample, Z — ► qq7 
decays were selected by tagging hadronic final states with isolated photon can- 
didates in the electromagnetic calorimeter. Combining the measured rates of 
Z — * qq7 decays with the total rate of hadronic Z decays permits the simultane- 
ous determination of the widths of the Z into up- and down- type quarks. The 
values obtained, with total errors, were 

T u = 300±1| MeV and T d = 381±^ MeV. 

The results are in good agreement with the Standard Model expectation. 
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1 Introduction 

Isolated hard photon production in hadronic Z decays provides information about the 
electroweak couplings of quarks and about QCD [1]. In e + e~ collisions at the Z res- 
onance, initial state radiation is suppressed and the production of isolated photons is 
mainly due to final state radiation (FSR) from the primary quark-antiquark pair. 

The aim of the analysis of hadronic events with FSR presented in this paper is to 
measure the electroweak couplings to up- and down-type quarks. The measurement is 
based on the photon's property of coupling to the electric charge. By selecting hadronic 
Z decays with FSR a subsample can be extracted which is enriched in up-type quarks 
since the strength of the coupling is proportional to the charge squared. Combining 
the measured rates of hadronic events with final state photon radiation with the total 
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rate of hadronic Z decays permits the simultaneous determination of the partial decay 
widths into up-type and down-type quarks. This approach is complementary to anal- 
yses of partial decay widths dealing with either heavy or light quark flavour tagging. 
This analysis determines the electroweak coupling of the up-type quarks u, c and the 
down-type quarks d, s, b for the flavour admixture present on the Z resonance. The 
LEP experiments have already published QCD measurements and electroweak quark 
coupling measurements based on final state photon production [2-10]. This analysis is 
based on a larger data sample and updates earlier OPAL measurements [8]. 

The Z decay into a quark-antiquark pair can be described by effective vector and 
axial-vector couplings, <7v,q an d #A,q, where q=u or d, representing up-type or down- 
type quarks. The Z decay widths depend on the sums of the squares of these couplings. 
Writing effective couplings c q as 

Cq = 4 (Sv,q + 0A,q) , (!) 

the total hadronic width of the Z, which is calculated to third order in a s [11], is given 
by: 



N r m%G 



a 2 a a 3 



r h adrons = "^Sif 1 + ~ + L4 ^ " 12 ' 8 ^ ( 2c u + 3d), (2) 



247rV2 V 71 71-2 tt 3 , 

where a s is the strong coupling constant, Gp is the Fermi constant, mz is the mass of 
the Z, and N c = 3 is the number of colours. The electromagnetic quark couplings also 
enter the hadronic partial width of the Z with a final state radiation (FSR) photon; 

N c myGF . . a /n 2 _ 2 N 
Thadrons+7 = F(y cut ) — (2q u c u + 3g d c d ), (3) 

where a is the electroweak coupling constant, q n and g d are the charges of up-type 
and down-type quarks, F(y cut ) is a correction factor determined from perturbative 
calculations, and y cnt is a jet resolution parameter, which determines whether a photon 
is merged into a nearby jet. These widths can be re-written in terms of decay widths 
of the up- and down-type quarks: 

r hadro „s = 2r u + 3r d; r hadrons+7 = — (l - ^ - L4 ^2 _ 12 .8a 3 /^) (8Fu + 3Td) - 

(4) 

The measurement of the partial hadronic width of the Z with an isolated FSR photon 
and the measurement of the total hadronic width of the Z depend on different linear 
combinations of c u and c d , or equivalently of T u and T d . In order to extract the partial 
widths or the couplings c Ujd , the factor F(y cat ) in equation 3 has to be calculated. 
Matrix element calculations at 0(aa s ) are available to perform these calculations as 
a function of the jet resolution parameter y cut [12]. One important aspect of these 
predictions is the regions of phase space in which they are valid. In particular, in- 
frared divergences due to soft photons and photons collinear with soft quarks must 
be avoided. For this reason, hard isolated photons are studied. A good agreement 
between hard isolated photon production and matrix element predictions has already 
been demonstrated by OPAL [13]. In the range of intermediate values of y cut , F(y cut ) 
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can be determined with a theoretical uncertainty of around 1% based on these predic- 
tions. The uncertainty is evaluated by varying a s and the phase space cut-off in the 
perturbative calculation as described in more detail in section 6. In contrast to the 
procedure followed in reference [9] no correction to account for the b-quark mass is in- 
cluded in the equations (2) and (3). The results of the present analysis were extracted 
at effective masses of the photon-jet system where the relative impact of the b-quark 
mass is much smaller than in the case of reference [9] . 

In this paper, we present a new measurement of the partial widths T u and T d . 
The full LEP1 data sample taken at the peak of the Z resonance was used in this 
analysis, amounting to more than 3 million hadronic Z decays. After describing the 
main characteristics of the OPAL detector in section 2, a description of the event 
selection follows in section 3. The background and efficiency corrections applied to 
the selected data are described in section 4. The systematic errors are discussed in 
section 5 before the results are presented in section 6. 

2 The OPAL Detector 

The OPAL detector operated at the LEP e + e~ collider at CERN. A detailed description 
of the detector can be found in reference [14]. For this study, the most important com- 
ponents of OPAL were the central detector and the barrel electromagnetic calorimeter. 
The central detector, measuring the momentum of charged particles, consisted of a sys- 
tem of cylindrical tracking chambers surrounded by a solenoidal coil which produced 
a uniform axial magnetic field of 0.435 T along the beam axis 1 . 

The electromagnetic calorimeters completely covered the azimuthal range for polar 
angles satisfying | cos#| < 0.98, providing excellent hermeticity. The barrel electromag- 
netic calorimeter covered the polar angle range | cos#| < 0.82. It consisted of 9440 lead 
glass blocks, each 24.6 radiation length deep, approximately pointing towards the inter- 
action region. Each block subtended an angular region of approximately 40 x 40 mrad 2 . 
Half of the block width corresponded to 1.9 Moliere radii. Deposits of energy in adjacent 
blocks were grouped together to form clusters of electromagnetic energy. The intrinsic 
energy resolution of <7e/E = 0.2% © 6.3%/\/E was substantially degraded (by a factor 
~ 2) by the presence of two radiation lengths of material in front of the lead glass. For 
the intermediate region, 0.72 < | cos#| < 0.82, the amount of material increased to up 
to eight radiation lengths. The two endcap calorimeters, each made of 1132 lead glass 
blocks, 22 radiation lengths deep, covered the region of 0.81 < | cos#| < 0.98. In this 
study the measurement of inclusive photon production is restricted to the barrel part 
of the detector. 

1 In the OPAL coordinate system, the x axis points towards the centre of the LEP ring, the y axis 
points approximately upwards and the z axis points in the direction of the electron beam. The polar 
angle 9 and the azimuthal angle <fi are defined with respect to the z and x-axes while r is the distance 
from the z-axis. 
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3 Event selection and reconstruction 



The first stage of the event selection was to identify hadronic Z decays. In this sam- 
ple, isolated, high energy, neutral clusters in the electromagnetic calorimeter were then 
considered as photon candidates. To reduce the dominant background from n° decays, 
the shower shapes were additionally required to be consistent with a shower produced 
by a single photon, as opposed, for example, to overlapping showers from the two pho- 
tons from a n° decay. Jets were then reconstructed in events with a photon candidate. 
The jet reconstruction implied an additional isolation requirement on the photons, as 
a function of the jet resolution parameter. 

Hadronic Z decays were selected from the data taken by the OPAL detector between 
1990 and 1995 at a centre-of-mass energy of 91.2 GeV. The data taken above and below 
the Z resonance during those years were not analysed. Using the criteria described in 
reference [15], A^z^Hadmns = 3 022 897 hadronic Z decay candidates were selected. The 
tracks and clusters were defined as described in reference [16]. 

Isolated electromagnetic clusters in the electromagnetic calorimeter were identified 
within the selected sample. The following criteria were used: 



1. Only clusters with no associated tracks were considered. 

2. The clusters had to be located in the barrel region of the electromagnetic calorime- 
ter defined by | cos Cluster | < 0.72. This requirement reduces the background from 
initial state radiation, which peaks at low polar angle. 

3. The energy of each cluster was required to be larger than 7 GeV. This reduces 
the background from neutral hadron decays and initial state radiation. It also 
enabled the use of a shower shape fit, for which the description of the Monte 
Carlo was not satisfactory below 7 GeV. 

4. The shape of the cluster was required to be consistent with a single photon. The 
number V b i oc ks of lead glass blocks contributing more than 0.02 GeV to the clus- 
ter energy was required to be less than 16. A cut at W < 30 mrad was set on 
the energy weighted first moment of the cluster, defined as: 




w= ^=1 , vv^ -r, v ; (5) 



blocks £meas _ ^2 + (Q. _ 0)2) 

E ^blocks Tpn 
1=1 

where E™ eas , (pi and 0, are the energy deposited in block i, its azimuthal and polar 
angles. and 9 are the and 9 positions of the cluster centroid. In addition, 
the shower-shape parameter C was required to be smaller than 5. The same 
definition of C as in [17] was used: 

c = -^—y {E \ ] } , (6) 
^blocks v w r ) 2 + (ar as ) 2 

where Ef md is the predicted energy in calorimeter block number %. This was taken 
from the best fit of the shower profile parametrisation to the observed energy 
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sharing between the calorimeter blocks assuming that the cluster was produced 
by an isolated photon. The reference profiles varied as a function of cos 9 because 
of the varying amount of material in front of the calorimeter, of red = 0.05 • Ef red 
is the error on Ef rcd where 0.05 is the accuracy of the Monte Carlo integration 
and a™ = 0.002 • E™ eas © 0.063^^^ is the estimated error on E™ eas [14]. 

5. The cluster candidates were required to be isolated. Clusters with tracks or other 
clusters of an energy sum greater than 0.5 GeV within a cone of half opening angle 
of 0.235 rad around the cluster were rejected. 

In total, 12 626 events with a photon candidate were selected from the original 
sample of hadronic Z decays. Each event in the sample was reconstructed into jets. 
Two jet finding schemes were used and compared: JADE E0 [18] and Durham [19]. 
Both are based on the resolution parameter 

M 2 

y^j = ( 7 ) 

where E vis is the total visible energy in the detector. In the Jade algorithm M i3 - is 
defined as the invariant mass 

Mfj JADE = 2^,(1 -cosc^-), (8) 

while in the Durham or kx algorithm it is defined as the minimum transverse momen- 
tum of one particle with respect to the other, 

M 2 Durham = 2 min(^ 2 , E|) ( 1 - COS « y ) , (9) 

where E iy Ej are the energies of the particles and is the angle between them. In 
each case, the jet finding was performed in two steps. In the first step, the jets were 
reconstructed excluding the photon candidate by iteratively combining particles until 
for all possible particle combinations the calculated yij was larger than a specified 
2/ cut . This was done for 12 values of y cnt in both algorithms. In the second step, y^i 
was calculated for the photon candidate and each jet i in the event. Each value of 
y^i was required to be greater than the y cut used during the first step. Events failing 
this requirement were rejected. The number of events n raw passing this last selection 
criterion depends on the chosen y cut value. The choice of values for the y cnt in each 
algorithm ensured that the photon candidate was sufficiently isolated from the jets and 
provided a scan of the phase space to allow a valid and more detailed comparison of 
the result with matrix element calculations. The selected sample size for different y cut 
values is shown in table 1 for the two jet finding schemes. 



4 Background and efficiency corrections 

The number of selected events was corrected for background contributions, detector 
response and selection efficiency Two categories of background were considered. The 
first was the contribution from misidentified neutral hadrons and their decay products. 
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The second was from genuine photons originating from initial state radiation. The 
sequence in which the data were corrected for background and efficiency reflects this 
difference in background sources. First, the background contribution from neutral 
hadrons was estimated by a fit to the shower shape distribution in data, using the 
difference in shower shape to distinguish genuine photons from hadron background. 
Second, a correction for detector effects was applied. In particular this corrected for 
the photon detection efficiency and compensated for the requirements on the shower 
shape and measurement uncertainties. It was calculated using a full simulation of the 
OPAL detector. Third, events with photons produced by ISR were subtracted. Their 
number was estimated with Monte Carlo simulations at the hadron level. Finally, the 
number of events remaining at the hadron level was corrected for hadronisation and 
selection efficiencies, so that it could be compared with results from matrix element 
calculations at the parton level. In this analysis, hadron and parton level were defined 
as in reference [20] . 

4.1 Neutral hadrons 

Neutral hadrons and their decay products can produce isolated clusters in the electro- 
magnetic calorimeter which pass the selection for photon candidates. This background 
was simulated using a sample of 8.4 million events generated with Jetset 7.4 [21], 
which were processed with the OPAL detector simulation program [22]. The hadronic 
background comes mainly from events with a 7r° decay, and has a small contribution 
from f] decays and K° L . There is also a contribution below 1% of the total hadronic back- 
ground from neutrons, oj and p. The simulated C distributions are shown in figure 1 for 
single photon clusters and for clusters produced by the hadronic background sources. 
By fitting the data distribution to a linear combination of the expected distribution of 
the shower-shape parameter C for the background and the measured distribution for 
single photons taken from a data sample of /i/i7 events, the fraction of single photon 
events was extracted [3]. In figure 2 the C distributions are plotted for isolated neutral 
clusters from /z/ry events in data and Monte Carlo and good agreement is observed. 
The fit is a binned maximum likelihood fit which takes into account the effects of data 
and Monte Carlo statistics [23]. For each y cut value considered, the fit was performed 
in eight cluster energy bins to extract the fraction of single photon events. The total 
fraction / 7 (y cu t) of single photon events for each y cut value was estimated by calculat- 
ing the weighted average of the results for the eight energy bins. The procedure was 
repeated for the two jet finding algorithms. Figure 3 shows an example of the fit results 
for the eight cluster energy bins. The jet finding scheme was in this case Jade E0, 
with y cut = 0.08. Table 2 shows the extracted fraction of single photon events in the 
selected samples for three of the y cut values in each jet finding scheme, spanning the 
phase space studied. The quoted 68.3% confidence level error was obtained by studying 
the variation of ln£ [23]. It combines the uncertainty due to limited data statistics 
(2526 fifi'j events, and between 1900 and 9000 candidates, see table 1) and due to 
limited Monte Carlo statistics (between 548 and 1695 clusters from neutral hadrons). 

The background estimate was found to be consistent with previous measurements 
performed by OPAL and L3. OPAL measured [24] that Jetset underestimates the 
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rate of isolated r] with an energy E v > 4.5 GeV by a factor 2.07 ±0.11, and the rate 
of isolated tc° in the energy range 4.5 < E n o < 22.8 GeV by a factor 1.99 ± 0.05. L3 
measured [25] a rate of isolated electromagnetic clusters which was 1.88 ± 0.08 times 
larger than predicted by Jetset. Applying the present method to our data sample 
over a range of cluster energies 4.5 < -E^ust < 46 GeV, we measured an excess of a 
factor 2.05 ± 0.03. The value is compatible with the previous measurements of OPAL 
and L3. 



4.2 Detector Effects 

After subtracting the events with an isolated neutral hadron from the sample of can- 
didates, the data were corrected for detector effects: 

n had(2/cut) = ^raw(ycut) X / 7 (y C ut) X C d< % cut ) (10) 

with n rsw (y cut ) as the number of photon events before the background correction, and 
f-y(Ucut) as the fitted single photon fraction (ISR+FSR) in the selected sample at detec- 
tor level, which has been derived as discussed in the previous section. The correction 
factor c det (y CVLt ), which takes into account detector inefficiencies, was calculated for 
each y cnt and each jet finding algorithm using the sample of 8.4 million fully simulated 
Jetset events. The same cuts on the photon energy, polar angle and isolation were 
applied at the hadron level. The correction factors c det (y cut ) are shown in table 2, again 
for three values of y cut in each jet finding scheme. The efficiency losses are mainly due 
to photon conversions in the two radiation lengths of material in front of the electro- 
magnetic calorimeter, e.g. the beam pipe, the central tracking system and the magnet 
coil. The photon conversion rate was calculated to be around 8%. 



4.3 Initial state radiation 

Initial state radiation photons are indistinguishable from photons from final state ra- 
diation. Although the energy and angular cuts on the electromagnetic clusters reduce 
the amount of initial state radiation, which is already small when working at the Z 
resonance, initial state photons can still make a contribution of several percent to the 
selected sample. This contribution was evaluated using the /C/C Monte Carlo [26] with- 
out the full detector simulation; the /C/C generator has the most accurate ISR modeling 
for e + e~collisions. The data which are corrected for detector effects are directly com- 
parable to the /C/C hadron level predictions. The number of FSR events at the hadron 
level is therefore given by 

nl!ny C nt)=nl ad (y cut )-b 1SR (y cnt ), (11) 

with b (|/cut) as the number of ISR events as calculated with the JCJC Monte Carlo 
program and normalized to the number Az^Hadrons of selected hadronic Z decays. The 
same cuts on the photon energy, polar angle and isolation were applied at the hadron 
level for both the Jetset and /C/C Monte Carlo samples. The numbers of isolated 
initial state photons b ISR (y cut ) to be subtracted are shown in table 2 for three values 
of y cut in each jet finding scheme. 
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4.4 Geometrical acceptance, isolation cone, and hadronisation 

In order to compare the measured final state photon rate with the matrix element 
predictions, the measurements must be corrected for the remaining acceptance cuts 
on the polar angle and the cone isolation, for hadronisation, and for fragmentation, 
which affect the number of jets in the event. A correction factor c pa,r (y cnt ) which 
takes into account the cuts mentioned above was calculated using Jetset 7.4. No 
further correction was made for the energy cut -E c iust > 7 GeV, to avoid introducing 
additional uncertainties. The matrix element predictions were calculated with the same 
requirement, E 1 > 7 GeV. The correction factors c par (y C ut) for three different y cut values 
in two jet finding algorithms are shown in table 2. The corrected numbers of isolated 
final state photon candidates per 1000 hadron events, 

^fsr = nH*(y cut ) x 1000 = ng?(Vc«t) x c^(y cut ) x 1000 , (12) 

* * Z— >Hadrons *Z— >Hadrons 

are shown in table 1. 



5 Systematic Effects 

The contributions to the systematic uncertainty are summarized in table 3. Relative 
errors are given for the value of y cnt in each scheme which minimizes the total uncer- 
tainty on the photon rate, and for the extreme y cut values considered. The statistical 
uncertainty quoted in table 3 combines the error from limited data and Monte Carlo 
statistics. It is calculated as the quadrature sum of the fit error (see section 4.1) and 
of the uncertainties of all Monte Carlo samples generated and used to simulate various 
other aspects of the analysis. The systematic studies from which the other uncertainties 
were obtained are described below. 



5.1 Background 

The estimate of the hadronic background which results from the fit of the shower-shape 
parameter C relies on the proper simulation of the shower for the background, as well as 
on the simulation of the yield of the various hadronic background sources. To quantify 
the systematic uncertainty originating from the description of the hadronic background, 
the fit procedure was repeated for all y cut values and the two jet finding algorithms after 
changing the Monte Carlo input distributions to which the data are fitted by varying the 
composition of the hadronic background. The contribution from n° and rj was doubled 
to account for the observed underestimation of isolated neutral pions and ^-mesons 
in the Monte Carlo [24]. The contribution from K° was increased by 18%, which 
was extrapolated from the observed differential cross-section in [27]. The observed 
differences from the results of the default procedure were not statistically significant, 
and showed no systematic trend, so no systematic uncertainty was assigned. Further 
systematic studies have been made by varying the bin size and the upper boundary of 
the C distribution. Again, no significant change in the resulting rate was observed. 
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The effect of the interference between initial state radiation and final state radiation 
was checked by comparing prompt photon samples from the /C/C Monte Carlo generated 
with and without ISR-FSR interference. The average difference of around 1 per mill in 
the rate of isolated prompt photons with an energy of larger than 7 GeV was assigned 
as systematic uncertainty on the corrected rate and is listed in table 3. 

5.2 Detector description 

Several systematic effects coming from the description of the detector were investigated: 

• The track parameters were smeared by 10% to estimate the effect of the recon- 
struction on the jet finding and the isolation criteria [28]. 

• The effect of the energy resolution of the electromagnetic calorimeter was eval- 
uated by smearing the energy of clusters by 20%. This number is chosen since 
the energy resolution of isolated photon clusters selected in Jetset and in a 
/z/ry reference data sample differ by up to 20%. The cluster energies were also 
shifted by 0.1-0.2 GeV. This shift reflects the difference in the mean values of the 
distributions of Cluster — E 1 from the simulation and of Cluster — Skin in data, 
where E kin is the kinematically reconstructed energy of the isolated photon in 
the reference sample. 

These effects would modify the correction factor c dct (y cut ). In all cases, the systematic 
uncertainty on the correction factor was found to be below the per mill level. 

An important contribution to the systematic error comes from the simulation of 
photon conversions. An error of 10% on the photon conversion rate [29] results in an 
error on the correction factor of about 0.8%. 

The jet reconstruction was also tested for systematic effects. The energy resolution 
was checked using a subsample of the candidate events with two jets and one isolated 
photon. The measured energy of each jet was compared to the energy of the same 
jet calculated by imposing energy/momentum conservation in the final state. The 
measured shift in the jet energy was used to estimate the systematic error on the 
photon rate. A variation between 0.1% and 0.6% in the rate of photons was found for 
the Jade E0 jet finding scheme as a function of y cut , and between 0.1% and 0.3% for 
the Durham scheme. 

The jet angle reconstruction was tested using the same sample of two-jet events. 
The hadronic part of each event in both Monte Carlo and data was boosted into its 
centre of mass using the photon momentum. In this system the two jets should be 
collinear. The angular shift between the Monte Carlo expectation and the data was 
taken to calculate a systematic uncertainty. This uncertainty varied between 0.1% and 
1.7% for both jet finders. Studies of the photon angular resolution showed that the 
effect on the rate was negligible. 

The combined systematic uncertainties obtained from these studies are given in 
table 3 as the contribution to the error from the determination of the correction factor 

C dct (ycut). 
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5.3 Fragmentation and hadronisation 

Three methods were used to evaluate the uncertainty arising from fragmentation and 
hadronisation. They affect the determination of the correction factor c par (y cu t) , which 
compensates for the geometrical acceptance, i.e. the extrapolation of the result to 
the total cos 6 range, and for the isolation cone around the photon candidate. An 
uncertainty is not trivial to assign since we largely rely on Monte Carlo studies. 



Parton shower model 

The Jetset predictions were compared with those of Herwig [30] and Ari- 
adne [31], which use different parton showering schemes. Herwig also uses a 
different fragmentation scheme. The relative deviations are larger towards the 
ends of the y cut range with a largest deviation of —9.0% (for Herwig and the 
Durham algorithm). The deviations are smaller in the intermediate range of 
y cut values for both jet finder schemes. The deviations stem from the different 
particle flow around the FSR photons in the three models. As a result, in the 
Jetset case more photons are rejected at the hadron level than in the Herwig 
or Ariadne cases due to the cone isolation requirement. This effect is smaller 
in the region of intermediate y cut values. 

First derivative prediction 

The photon rates measured with different y cut values are correlated. The first 
derivative of the photon rate is almost statistically independent for the different 
y cut values. Comparing the first derivative distribution for different parton shower 
models provides an additional systematic cross-check of the efficiency correction. 
The first derivative is defined as 

D FSR {y cnt ) = -tt — T[N FS R{ycnt-S/2)-N FSR {y cut + 5/2)] 

^Z-^Hadrons^ 

1 jyi^rcjcct _ jyrcjcct-^ij 

Z^Hadrons^ 

1 1 

where Afsr is the number of final state photon candidates. j\r* ^reject j g 
number of events that have a jet multiplicity i at y cnt — S/2 but that are rejected at 
y C ut + <V2, whereas jV reject ~^ is the number of events that are rejected at y cut — 5/2 
but are retained at y cnt + 5/2 with a jet multiplicity i. The gains and losses for 
each i-jet event class are statistically independent, so that the quantities, 

xj-j— >reject ly-reject— ►« 

c «^reject _ parton ^ c reject->i _ parton q^n 

AT-i— >reject ' 1 » 7-rejcct^i \ / 

hadron hadron 

are also almost statistically independent. The discrepancy between these quanti- 
ties calculated for Herwig and Jetset is largest for small y cut values. Relative 
to the Jetset value, it amounts to up to -18.5% for the Jade E0 scheme and 
up to -24.1% in the Durham scheme. For y cut > 0.02, the predictions agree to 
within 6% for the Jade E0 scheme. For the Durham scheme, they agree within 
10% for y cut > 0.008. The deviations are also due to the cone isolation criterion. 



N, 
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• Isolation cone 

The analysis was repeated with different isolation cones, varying between 0.175 rad 
and 0.335 rad and allowing for up to 0.5 GeV of additional energy within the iso- 
lation cone. The largest deviations are observed for a cone of 0.175 rad half angle. 
For the Jade E0 jet finder, they are below 5% for all values of y cut and have a peak 
value of 4.9% for y cut = 0.2. The deviations calculated in the Durham scheme 
are at the same level, with the highest value of 4.8% at the largest y cut = 0.1. 

Because the previous considerations are not independent, the estimates which give 
the largest positive and the largest negative contributions to the systematic uncertainty 
for each y cnt have been assigned as systematic errors on the photon rate. They are listed 
in table 3 as the error contribution from the determination of c par (|/ cut ) and they are 
typically largest for small y cut values, smaller for intermediate y cut values, and increase 
again for large y cut values. This indicates that the correction c par (y cut ) and therefore the 
final result is more sensitive to details of the energy flow around the photon candidate 
for very large and very small values of y cut , i.e. near the phase space boundaries. 



Figure 4 shows the measured rate of FSR photons normalised to 1000 hadronic Z 
decays as a function of the y CVLt value for the two jet finding algorithms. The results are 
compared to bands of matrix element Monte Carlo predictions for 0.15 < ai 1 ^ < 0.22 
and 0.0005 < y 7 < 0.001. is the first order a s value and y 1 is the cut-off in the 
phase space of the photon emission which is used in the matrix element calculations to 
remove the singularity of the 1/57,1/57 poles [4,8, 12,13]. The lowest two values of y cut in 
both jet finder schemes are not considered because of the large contribution from events 
with more than three jets, which are not described by the matrix element calculations. 
For both jet finding algorithms, the measurement errors in the intermediate y cut range 
are smaller than in the regions of small and large y cut values. At small and large y cnt 
values, the result is particularly sensitive to the isolation cone requirement, causing a 
large efficiency loss and a large systematic error. 

The decay widths T u and r d can be calculated from the measured FSR photon 
rate. The value y cut = 0.08 for the Jade E0 jet finding algorithm and y cut = 0.016 for 
Durham were chosen for the final results since they are the values for which the error on 
the FSR photon rate is the smallest. Taking a value of ai 1 ^ = 0.18±0.03 [32], correction 
factors F(y cut ) of 9.60l^ (y cut = 0.016, Durham scheme) and 6.11±g;g| (y C ut = 0.08, 
Jade E0 scheme) were calculated to be used as input to equation 3. With values of 
a s = 0.1172±0.002 and r hadrons = 1744.4±2.0 MeV [33] entering equation 2, we obtain 



6 Results 



8r u + 3r d = 3.55 ±0.08 



+0.07+0.03 
-0.07-0.01 



GeV 



(15) 



for the JADE E0 jet finding algorithm [y, 



cut — 



0.08) and 



8r u + 3r d = 3.50 ± 0.09±8;i2±8:oe GeV 



(16) 
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for Durham (y cut = 0.016). The first error is due to statistics, the second error is due to 
systematics and the third error comes from the uncertainty in ai 1 ^ and the theoretical 
cut off y T The results derived with the Jade and Durham schemes are consistent and 
the result with the smaller uncertainty (Jade scheme) is considered the main result. 

Figure 5 shows the correlation between the Z decay widths in up-type and down- 
type quarks obtained from this measurement. Also shown are the Standard Model 
prediction and the correlation obtained from the measurement of r hadrons = 2r u + 31^. 
Combining the information on Thadrons+7 and r hadrons, the decay widths for up-type 
and down-type quarks can be extracted for y cut = 0.08 in the Jade scheme: 

r u = 300 ± 13t^ MeV and T d = 381 ± 9^ MeV. (17) 

Note that in this measurement T u and T d are 100% anti-correlated. The Standard 
Model predictions [33] for these widths are: 

T u = 300.3 ±0.14 MeV and T d = 380.7 ± 0.15 MeV (18) 

in good agreement with the measurement. The result is also in good agreement with 
previous measurements from DELPHI [10], L3 [9], and OPAL [8,34] and is significantly 
more precise. The analysis in [34] used high-momentum stable particles as a tag, while 
the other three were based on the determination of the FSR rate in hadronic Z decays. 



7 Conclusion 



Using the entire OPAL LEP1 on-peak Z hadronic decay sample a measurement of the 
partial decay widths of the Z into up- and down-type quarks was performed. The 
measurement was based on 12 626 observed candidate hadronic Z decays with final 
state photon emission. The partial width Thadrons+7 for such decays was calculated 
as a function of the partial widths of the Z into up- and down-type quarks with the 
aid of a matrix element Monte Carlo program of 0(aa s ). The calculated value for 
Thadrons+7 was combined with the PDG average of the total hadronic Z width r hadrons 
to determine T u = 300l}g MeV and T d = 38lt\l MeV - The result is consistent with the 
Standard Model expectation and with previous measurements [8-10,34]. Its precision 
is significantly improved compared to previous measurements using a similar method of 
selecting hadronic Z decays with highly isolated and energetic photons due to increased 
data and Monte Carlo statistics. With respect to the previous OPAL measurement [8] 
a different approach was chosen to estimate the background rate of neutral hadrons 
being misidentified as genuine photons. The previous measurement invokes isospin 
symmetry to estimate the background rate of neutral hadrons and measures the rate 
of isolated charged pions to extract the rate of isolated neutral pions. The present 
analysis utilises the different shapes of the distribution of the shower-shape parameter 
C for neutral hadrons and genuine photons. The larger data and Monte Carlo statistics 
allowed us to choose this method which exhibits a smaller systematic uncertainty. 
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Jade EO 


Durham (k t ) 


Vcut 


^raw 


-Rfsr 


Vcut 


^raw 


-Rfsr 


0.005 


8503 


4.22 ±0.06 


0.002 


8936 


4.55 ±0.07 


0.01 


7504 


3.55 ±0.06 


0.004 


8017 


3.96 ±0.06 


0.02 


6351 


2.93 ±0.05 


0.006 


7331 


3.53 ±0.06 


0.04 


5039 


2.31 ±0.05 


0.008 


6876 


3.25 ±0.06 


0.06 


4140 


1.89 ±0.04 


0.01 


6366 


3.04 ±0.05 


0.08 


3496 


1.54 ±0.03 


0.012 


5996 


2.83 ±0.05 


0.1 


2904 


1.29 ±0.03 


0.014 


5613 


2.53 ±0.08 


0.12 


2513 


1.07 ±0.04 


0.016 


5331 


2.39 ±0.06 


0.14 


2261 


0.91 ±0.04 


0.02 


4890 


2.15 ±0.06 


0.16 


2094 


0.81 ±0.05 


0.04 


3743 


1.54 ±0.05 


0.18 


1991 


0.72 ±0.05 


0.06 


3288 


1.32 ±0.06 


0.2 


1933 


0.68 ±0.05 


0.1 


2913 


1.23 ±0.06 



Table 1: Total number of candidate events n raw and corrected number of selected 
candidate events in 1000 hadron events -Rfsr m the two jet finder schemes for twelve 
different values of y cut . The error represents the combined uncertainty from data and 
Monte Carlo statistics. 





Jade E0 scheme 


Vcut 


0.005 


0.08 


0.20 


f-y(y 'cut) 


0.851 ±0.011 


0.865 ±0.016 


0.670 ± 0.042 


c dct {y C ut) 


1.142 ±0.005 


1.152 ±0.008 


1.073 ±0.011 


& 1Si %cut) 


354 ±8 


183 ±6 


83 ±4 


c par (2/cut) 


1.612 ±0.008 


1.413 ±0.010 


1.573 ±0.018 




Durham (kr) scheme 


Vcut 


0.002 


0.016 


0.1 


ff (ycut) 


0.839 ±0.012 


0.848 ±0.019 


0.719 ±0.035 


c dct (y cu t) 


1.146 ±0.004 


1.137 ±0.006 


1.062 ±0.008 


b^(y cu t) 


380 ±8 


222 ±6 


109 ±5 


c par (t/cut) 


1.673 ±0.008 


1.471 ±0.008 


1.754 ±0.017 



Table 2: The single-photon fraction / 7 (y cu t), the correction factor c det (y CVLt ) which takes 
into account detector inefficiency, the expected number b 1SR (y cut ) of background events 
due to initial state radiation in the sample of selected candidates, and the efficiency 
correction c par (y cut ). For each jet finding scheme, values are shown for the y cut which 
gives the minimum total error, and for the two extreme y cut values. The errors reflect 
the available Monte Carlo statistics and in the case of the single-photon fraction the 
combined data and Monte Carlo statistics. 
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Relative 


error in % for 










Jade EO scheme 


Durham (kx) 


scheme 


Vcut 


0.005 


0.08 


0.20 


0.002 


0.016 


0.1 


Statistics 


±1.5 


±2.2 


±6.8 


±1.6 


±2.4 


±5.3 


(candidates, /z/ry, MC) 














Systematics 














ISR background 


±0.1 


±0.1 


±0.1 


±0.1 


±0.1 


±0.1 


C d6t (l/cut) 


±1.8 


±1.2 


±0.8 


±1.8 


±1.5 


±0.8 


C par (2/cut) 


±1.5 


±1.6 


±4.9 


±1.8 


±2.5 


±4.8 




-18.5 


-1.6 


-3.7 


-24.1 


3.6 


-9.0 


Sum 


±2.7 


±3.0 


±8.4 


±3.0 


±3.8 


±7.2 




-18.6 


-3.0 


-7.8 


-24.2 


-4.6 


-10 



Table 3: Summary of error contributions for the cases with the smallest total un- 
certainties in each of the two jet finder schemes. The error contributions for the two 
extreme values are also shown. The estimation of the various contributions is discussed 
in section 5. 
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Shower-shape parameter C 








1 2 3 4 5 
Shower-shape parameter C 



Figure 1: The left plot shows distributions of the shape variable C for single photon 
and neutral hadron clusters as obtained with the OPAL detector simulation program. 
The right plot shows the simulated C distributions for various hadronic background 
sources in the case of the Jade scheme for a value of y cnt = 0.08. 
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Figure 2: The distributions of the shape variable C for photon clusters in /i + /i~7 events 
for data and Monte Carlo classified in eight bins according to their measured energy. 
The Monte Carlo distributions are normalised to the data statistics. 
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Figure 3: The data distributions of the shower-shape parameter C in eight bins accord- 
ing to the measured cluster energy and the corresponding fit result in each bin. The 
plots show the results in the Jade E0 scheme for y cnt = 0.08. The filled histograms 
represent the distributions of the neutral hadron background scaled to the fractions 
resulting from the fit. The open area bordered by the line represents the single photon 
contribution on top of the hadronic background. Note that for E c \ > 40 GeV (last 
plot) the best fit only yielded a very small contribution from hadronic background. 
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Figure 4: The total number of events with FSR per 1000 hadron events in two jet finder 
schemes and for twelve values of y cut . The prediction of the 0(aai^) matrix element 
Monte Carlo is also shown. The band width of the theoretical prediction is a result of 
a variation of ag and of the phase space cut-off in the perturbative calculations. 
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Figure 5: Correlation between the Z decay widths to up-type and down- type quarks. 
The band width of Thadrons+7 = 8r u + 31^ reflects one standard deviation as obtained 
from this measurement. The narrow band shows the correlation within one standard 
deviation as obtained from the world average of r hadrons = 2r u + 3r d . The point on 
the narrow band represents the Standard Model expectation. 
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